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Abstract 

We perform an analytical study of the correspondence between a 
classical oscillator with frequency perturbed by a coloured noise and 
the one-dimensional Anderson-type model with correlated diagonal 
disorder. It is rigorously shown that localisation of electronic states in 
the quantum model corresponds to exponential divergence of nearby 
trajectories of the classical random oscillator. We discuss the relation 
between the localisation length for the quantum model and the rate 
of energy growth for the stochastic oscillator. Finally, we examine the 
problem of electron transmission through a finite disordered barrier 
by considering the evolution of the classical oscillator. 

Pacs numbers: 05.40.-a, 71.23.An, 72.15.Rn 

1 Introduction 

This work serves the goal of establishing some quantitative links between 
two seemingly unrelated fields: quantum disordered models on the one hand 
and classical stochastic systems on the other. More precisely, we analyse the 
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relations existing between a classical oscillator with frequency perturbed by a 
feeble noise and the one-dimensional (ID) Anderson-type model with a weak 
diagonal disorder. Our main interest is in correlated random potentials and, 
correspondingly, in coloured noise for the stochastic oscillator. 

Recently, the role of correlations in random potentials of quantum mod- 
els has been the object of intense scrutiny. In particular, it was shown that 
specific long-range correlations in potentials may lead to the emergence of 
a continuum of extended states even in ID lattices (see, e.g., [I], and ref- 
erences therein). In this paper we show that the phenomenon of Anderson 
localisation has its counterpart in the energetic instability of a random oscil- 
lator. Specifically, the mobility edge generated in the ID quantum models by 
long-range correlations is equivalent to the suppression of the energy growth 
of the stochastic oscillator due to temporal correlations of the frequency noise. 

We use the correspondence between stochastic oscillators and disordered 
solid state models in order to study the transmission properties of finite 
lattices by making use of the dynamical analysis of an oscillator with noisy 
frequency. This approach allows us to put in a new perspective the problem 
of electronic transport in disordered lattices and also to gain new insight on 
the dynamics of random oscillators. 

This paper is organised as follows. In the following section we define the 
models that constitute the object of our study, and we make some general 
considerations on their analogies. In Sec. ^|we rigorously analyse the relation 
between the localisation of electronic states for the Anderson model and the 
orbit instability of a random oscillator. In Sec. |] we discuss how correlations 
of the frequency noise can suppress the energy growth of the stochastic os- 
cillator. The analogy between a random oscillator and a disordered chain 
is then used in Sec. |5] to study the electronic transmission through a finite 
disordered lattice. In Sec. [5] we also discuss the relation between energetic 
instability and orbit divergence for a random oscillator. Finally, Sec. § is 
devoted to the summarising conclusions. 

2 Definition of the models 

The Anderson model is defined by the discrete stationary Schrodinger equa- 
tion 

tpn+l + 4>n-l + £n?Pn = Elp n (1) 
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where ijj n is the amplitude of the wave function at the nth site of the lattice 
and disorder is introduced via the site energies e n which are assumed to 
be random correlated variables. We do not restrict our considerations to a 
specific distribution for the random potential e n ; we only suppose that it has 
zero average (e n ) = and that the binary correlator (e n e n+ k) is a known 
function of the index k. We also assume that the correlator (e n e n+ k) does 
not depend on n and that it is a decreasing function of k. In other words, we 
make the physically sensible assumptions that the random succession {e n } is 
stationary, and that correlations decay with increasing distance. We restrict 
our analysis to the case of weak disorder, defined by the condition 

(el) « 1. 

In the preceding expressions the symbol (...) stands for the average over a 
single disorder realization defined by the limit 
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n=l 

we assume that this average is equivalent to the average over disorder real- 
izations (ensemble average) for the succession {e n }. 

It is known that the model (|I|) can be put into correspondence with the 
kicked oscillator defined by the Hamiltonian 

# = ^ (y + y ) + y ( A n 8(t - nT)) , (2) 

which represents an oscillator whose momentum undergoes instantaneous 
variations of random intensity A n at regular time intervals. The connection 
between the models (|1|) and (0) has been discussed before (see, e.g., [||]). 
Basically, the correspondence consists in the fact that, by integrating the 
Hamilton equations of motion of the oscillator (|2|) over the period between 
two successive kicks one gets the map 

x n+ i = x n cos(^T) + (p n - A n x n ) sin(^T) , . 

Pn+i = -x n sin(u;T) + (p n - A n x n ) cos(cjT) 

where x n and p n stand for the position and momentum of the oscillator 
immediately before the nth kick. This map is equivalent to to the Schrodinger 



3 



equation ([]]) which defines the Anderson model. Indeed, by eliminating the 
momentum from Eqs. (Q), one gets the relation 

x n+ i + x n _i + A n sm{u)T)x n = 2x n cos{uiT) 

which coincides with the Schrodinger equation (|l|) provided that the position 
x n of the oscillator at time t = nT is identified with the electron amplitude 
ip n at the nth site and that the parameters of the kicked oscillator are related 
to those of the Anderson model by the identities 

e n = A n sm(uT) and E = 2cos(uT). (4) 

The formal correspondence between the quantum model (|l|) and the 
kicked oscillator (Q) raises the question of whether a similar analogy can link 
the Anderson model to a random oscillator whose frequency is perturbed by 
a continuous noise instead that by a succession of discontinuous and singular 
kicks as in model (^|). In other words, one is led to infer the existence of close 
ties between the quantum model (P and a stochastic oscillator defined by 
the Hamiltonian 

"=«(t + t) + t«'> < 5 > 

where £(t) is a continuous and stationary noise. Notice that these require- 
ments on set the random oscillator fl5|) and the kicked oscillator (0) in 
two different categories within the vast family of stochastic oscillators, since 
the succession of kicks in the model (|2]) is a non-stationary and strongly 
discontinuous random process. Consequently, the connection between the 
Anderson model ([[[) and the kicked oscillator (|2|) does not prove at all the 
equivalence of models (|l|) and @ but only constitute a hint that such a link 
may exist. 

In our analysis of stochastic oscillators, we will focus on the Hamiltonians 
represented by Eq. (Q), completing the definition of the model by further 
assuming that the noise £(t) has zero average and that its binary correlator 
is a known function 

(£(*)> = and + t))= X (t). (6) 

In Eq. @ the symbol (...) is used for the time average 

(/(f)) = Hm 1 f T ° fit) dt, 

T ->OO In Jo 
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which is assumed to coincide with the ensemble average for the process £(t). 
Notice that we do not restrict our consideration to the case of white noise, 
but we are instead interested in the general case of coloured noise. Finally, 
we require that the noise £(£) be weak; in other words, we assume that the 
fluctuations of the frequency around its average value are small. 

Below we show that oscillators of the kind (|5|), with the above-mentioned 
noise features, are equivalent to the Anderson model (|I|) if two further con- 
ditions are met. First, the correlation function has to be of the form 



/ A2\ +00 

X{r) = hf £ ak)5(r-kT), (7) 

k=— oo 

where the symbol ((k) stands for the normalised binary correlators 

of the random variables A n specified by the second condition. Our second 
requirement is that the unperturbed frequency uj of the oscillator and the 
parameters A n must be related to the parameters E and e n of the Anderson 
model through the identities (f|). 

Notice that the links established by these two conditions associate key 
features of the noise £(t) to the corresponding properties of the random po- 
tential e n . Indeed, once the random variables e n and A n are connected by 
the relation (^), the correlators (§) become identical to the normalised cor- 
relators of the potential e n . Therefore the spatial correlations of the disorder 
in the Anderson model are mirrored by temporal correlations for the noise 
£(£). In the special case in which the disorder in the Anderson model is un- 
correlated (i.e., (e n+ k^n) — for k ^ 0), the noise for the random oscillator 
is white (i.e., (£(£)£(£ + t)) oc S(t)). One can also observe that the case of 
weak disorder in the Anderson model corresponds to that of weak noise for 
the random oscillator, since the condition (e^) <C 1 entails the consequence 
that (Al) <C 1 (except that at the band edge, i.e., for ujT — > 0, which is a 
special case where anomalies are expected to arise and will not be considered 
here) . 

Obviously, we must endow with a well-defined meaning the notion of 
"equivalence" used above to describe the connection between the Anderson 
model (|^) and the random oscillator defined by Eqs. (|3p and (0). We speak 
of equivalence of the two models in the sense that the time evolution of the 
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orbits of the random oscillator closely mirrors the spatial variation of the 
electronic states on the lattice. More precisely, the exponential divergence 
rate of nearby orbits turns out to be equal to the inverse localisation length 
of the Anderson model. 

The correspondence between the random oscillator ([5]) and the Anderson 
model (fl|) is to some extent surprising since the former is a classical system 
and is continuous in time whereas the latter model is quantum and discrete 
in space. It is therefore particularly interesting to notice how close the two 
systems turn out to be. To sum up, one of the main results of this paper is 
that the Anderson model with weak correlated disorder has a close analogue 
in a random oscillator with frequency perturbed by a coloured noise. This 
equivalence generalises the result established in Ref. j3j where the Anderson 
model with uncorrelated disorder was linked to a random oscillator of the 
kind (^) with white noise. 

3 The Lyapunov exponent 

In the previous section we have described the analogy between the Anderson 
model ([I]) and the random oscillator (|5|) as being based on the correspon- 
dence between the electronic wave-function of the former model and the space 
orbits of the latter system. To prove this analogy, we will compute the diver- 
gence rate of nearby trajectories of the random oscillator, i.e., its Lyapunov 
exponent and we will show that, when the conditions and (|7]) are met, 
the Lyapunov exponent coincides with the inverse localisation length in the 
Anderson model. We define the Lyapunov exponent through the formula 

,. ,. 11 f T \ x(t + 6) , 
A = lim lim— - / In y , > dt. 9 

To compute this expression it is convenient to introduce the polar coor- 
dinates defined through the standard relations x = rsin#, p = rcosO. This 
allows one to cast Eq. (|^) in the form 

1 [To -f 

A = lim — / - dt. 

T -Kx> T Jo r 

To proceed further, we consider the dynamical equations for the random 
oscillator in polar coordinates 

= uj + £(t) sin 2 0, (10) 
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r = -ir£(*)sin20; (11) 

Using the radial Eq. flTT|), the expression for the Lyapunov exponent can be 
finally put into the form 

A = - lim -L f T °at)sm(2d(t)) dt = -i(£(t) sin (26(t))). (12) 

T ->oo ZIq JO 2 

The problem of computing the Lyapunov exponent @ is thus reduced to 
that of calculating the noise-angle correlator that appears in Eq. (fL2|). This 
can be done in the following way, which is the extension to the continuum case 
of the procedure adopted in [|IJ for the discrete case. First, one introduces 
the noise-angle correlator defined by the relation 

7 (r) = <£(*)exp(2i0 (t + r))). 

Starting from this definition, in the limit e — > one has 

7(t + e) = (£(*) exp (z20 (t + r)) (l + 2i# (t + r)e)) + o(e). 

With use of the dynamical equation (pTOf ) one can further write 

7 (r + e) = 7 (r) (1 + 2icue) 

+2ie(f (t)f (t + t) exp (2i0 (t + r)) sin 2 (t + r)) + o(e). 

In the limit of weak noise, one can factorise the correlator that appears in 
the right hand side of the preceding equation and take the average over the 
angular variable using a flat distribution for 6. Indeed, when £(£) — > 0, 
Eq. (|T0|) implies that 9 ~ uj so that, after a conveniently long time, one 
can expect the angular variable to take values uniformly distributed in the 
interval [0, 2tt]. As a consequence the noise-angle correlator must obey the 
relation 

7 (r + e) = 7 (r) (1 + 2iue) - % - X (r)t + o(e), (13) 

where x{ T ) is the correlation function (or noise-noise correlator) defined by 
Eq. (^]). On the other hand, a simple application of calculus rules leads to 

d'y 

7 (r + e) = 7 (r) + ^(r) e + o(e). (14) 

Comparing Eqs. ( |i~3f ) and fll4l) , one obtains the differential equation 

d'y i . N 

— (r) = 2^ 7 (r) - -x(r), 
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whose solution (with the boundary condition lim r ^„ 00 j(t) = 0) gives the 
noise-angle correlator 



7 (t) = --[ x{s)e 2 ^ r - s) ds. 
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Using this result the Lyapunov exponent fljjD can be finally written as 

A=a/ (e(*)e(* + r))cos(2o;r)dT (15) 

o J — oo 

which implies that the Lyapunov exponent for the stochastic oscillator ([5]) 
is proportional to the Fourier transform x(2o>) of the correlation function at 
twice the frequency of the unperturbed oscillator. 

We are interested in the particular case in which the correlation function 
of the noise takes the specific form (^), because we want to prove that 
in that case the Lyapunov exponent fllBD coincides with the inverse localisa- 
tion length of the Anderson model ([I]). The substitution of the correlation 
function (0) in the general expression ( |l5l) gives 



+oo 



2£C(*0 cos(2u/77c) 



k=l 



Taking also into account the relations (f4|) between the parameters of the 
systems ([!]) and fl5|), one can finally write the Lyapunov exponent for the 
random oscillator as 

1 (e 2 ) +oc 
\ = -—^f—y{uT) with V (uT) = l + 2j2ak)cos(2ujTk). 
T 8 sm z (uT) ^ 

(16) 

This expression coincides with the one given in U for the localisation 
length in the Anderson model with correlated disorder. The inverse locali- 
sation length is given by the product of two factors, namely the Lyapunov 
exponent for the uncorrelated disorder case and the function (p(uT), which 
describes the effect of disorder correlations (and which reduces to unity when 
correlations are absent). Formula (|IB|) thus confirms the equivalence of the 
quantum Anderson model ([!]) with the classical oscillator (|5|) which had been 
inferred in Sec. ^| by the existence of a third system -the kicked oscillator (|2D~ 
which was somehow contiguous to both model ([!]) and (^). To sum up, for- 



mula (16) allows one to conclude that the Anderson model with correlated 
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disorder has a classical counterpart represented by a stochastic oscillator 
with frequency perturbed by a coloured noise. This conclusion generalises 
the equivalence established in j|] between the Anderson model with uncor- 
related disorder and an oscillator with frequency perturbed by a white noise. 

A remark is in order here: expression fll6|) for the inverse localisation 
length of model (|l|) is correct for all energy values inside the unperturbed 
band except that at the band centre, i.e., for uT = tt/2 where a anomaly 
arises and special methods are required for the analytical investigation (see, 
e.g., H). This anomaly is a resonance effect inherent in the discrete nature 
of the model (|1|) and cannot therefore be reproduced by the continuos sys- 
tem (H). Other anomalies appear in the Anderson model for the "rational" 
values of the energy (i.e., when uoT = 7ip/2q with p and q integer numbers), 
but they are effects of order higher than the second || and need therefore 
not be considered here. In conclusion, apart from the exceptional case of the 
band centre, the dynamical features of the models ([!]) and (|5[) do not differ 
to the second order of perturbation theory. 

The equivalence of the models (|IJ) and (|5]) can be examined also from 
a different point of view: that of the correspondence between discrete and 
continuous solid state systems. Indeed, the dynamical equation for the oscil- 
lator © 

x + u£(t)x = — to 2 x (17) 
coincides, mutatis mutandis, with the stationary Schrodinger equation 

- if;" + k£(x)i) = k 2 i) (18) 

which describes the motion of a quantum particle of energy E = k 2 in a 
random potential v(x) = k£(x). Actually, expression fll5|) for the inverse 
localisation length has long been known to solid state physicists (see, e.g., 
Ref. H) as the high-energy limit of the Lyapunov exponent for the continu- 
ous model fll8|) . Thus, the deduction of the inverse localisation length fll6|) 
for the Anderson model from expression (|l^) may be interpreted as the proof 
that the continuous model flT8"| ) can be put into one-to-one correspondence 
with the discrete lattice (0) if, and only if, the correlation function of the 
random potential has the specific form (|?p. (Obviously, the transposition of 
results from one model to the other requires a proper change of the corre- 
sponding parameters with relations like (|J); as a consequence of this swap, 
the mathematical correspondence of the two models does not imply an ex- 
act physical equivalence. Models ([!]) and (18), for instance, have different 
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unperturbed energy spectra, defined by the respective dispersion relations 
E = 2 cos k and E = k 2 ). 

4 "Mobility edge" for a stochastic oscillator 

In Ref. the authors used formula ([16]) to investigate the problem of mobil- 
ity edge for the Anderson model (|IJ). They showed that long-range correla- 
tions in the disorder can generate a continuum of extended electronic states 
and they found a way to construct sequences {e n } of site energies giving 
rise to a Lyapunov exponent with pre-defined dependence on the energy. In 
particular, using this recipe they were able to construct site potentials that 
generate a mobility edge even for the ID lattice ([I]). 

Here, we show how it is possible to solve the analogous problem for the 
random oscillator (|5|) taking formula (p~5|) as a starting point. More precisely, 
we will show how to define a continuous noise £(t) such that the correspond- 
ing Lyapunov exponent A (a;) has a pre-definite dependence on the frequency 
uo. Since the Lyapunov exponent determines the asymptotic behaviour of the 
oscillator energy (we discuss this point more in detail in the next section), 
shaping the function A(o>) through noise control enables one to determine 
the energetic behaviour of the oscillator. In particular, if the noise £(£) has 
the appropriate time correlations, the corresponding Lyapunov exponent can 
sharply drop from positive values to zero when the unperturbed frequency 
uo crosses a threshold value. In physical terms that means that the ener- 
getic growth of the oscillator is suppressed when the frequency reaches a 
critical value. The existence of a frequency threshold determining whether 
the oscillator is energetically stable or not is the physical counterpart of a 
mobility edge, which divides extended states from localised ones in the An- 
derson model. Thus, in spite of the current wisdom that frequency noise 
produces energetic instability (see, e.g., J?| and references therein), it turns 
out that time correlations of the noise may lead to a suppression of the en- 
ergy growth. This conclusion follows directly from the known formula ([To]), 
but, to the best of our knowledge, this implication has not been discussed 
before in the literature. 

To construct a noise £ (t) that gives rise to a defined Lyapunov exponent 
A(a>), the starting point is the correlation function x( T ) that can be easily 
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obtained by inverting formula ( |H>|) 

X (r) = - / XHe 2 ^ du. 



Once the correlation function x( r ) is known, we can obtain a stochastic 
process £ (t) satisfying the conditions ([]) by means of the convolution product 

/+oo 
(3(s) V (s + t)ds, (19) 
-oo 

where the function f3(t) is related to the Fourier transform x{ u ) °f t ne noise 
correlation function through the formula 

and 77 (i) is any stochastic process such that 

<77(t)> = and (ri(t)r}(t')) = 8(t - £')■ (20) 

Formula ([19]) defines the family of noises corresponding to a specific form A (a;) 
of the frequency-dependent Lyapunov exponent and constitutes the solution 
to the "inverse problem" (i.e., determination of a noise that generates a 
pre-defined Lyapunov exponent). 

As an example, we can consider the Lyapunov exponent 

AM = { 1 if < 1/2 (21) 
[ otherwise ' 

whose frequency dependence implies that the random oscillator undergoes a 
sharp transition for \u\ = 1/2, passing from an energetically stable condition 
to an unstable one. Following the described procedure it is easy to see that 
the Lyapunov exponent fl2T|) is generated by a noise of the form 

"+00 



V» f +oc sin Is) . \ , 

IT J -00 S 



with r](t) being any random process with the statistical properties fl20|). 

At this point, it is opportune to stress that the mathematical identity of 
Eqs. (^) and (|TH) implies that all features of the random oscillator (|5|) are 
shared by the solid state model (0). Therefore the mathematical results of 
this Section not only imply that noise correlations can make the random oscil- 
lator (^|) stable; they also represent a recipe to construct a random potential 
generating a mobility edge for the model (M 
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5 Transmission through a disordered barrier 



We are now in the position to see how the analogy between the quantum 
model (HD and the random oscillator @ can be used not only to compute 
the localisation length in the Anderson model but also to deal with problems 
both more challenging and of greater physical interest, such as the study of 
the transmission properties of a disordered barrier. In this Section we show 
how the random oscillator formalism allows us to tackle this problem and 
how it is possible to obtain expressions for the transmission coefficient as a 
function both of the sample length and of the inverse localisation length fll5|). 

More specifically, let us consider the case of a ID disordered lattice of L 
sites sandwiched between two semi-infinite perfect leads. Mathematically the 
problem is defined by the Schrodinger equation ([I]), where the site energies 
e n are now equal to zero for n < 1 and n > L, while for 1 < n < L they 
are assumed to be correlated random variables. In || it was shown that the 
transmission coefficient Tj, through the L-sites segment can be expressed in 
terms of the classical map (||) as 

4 

Tl = 2 + r 2 (L)+r 2 (L)' (22) 

where ri(L) and r2(L) represent the radii at the Lth step of the map tra- 
jectories starting from the phase-space points Pi = (xo = l,Po = 0) and 
Pi = (^o = 0, po — 1) respectively. An analogous formula was given in || for 
continuous models like the one defined by Eq. (|18|). 



Formula (|22|) constitutes the bridge that makes possible to link the trans- 
mission properties of a disordered barrier to the time evolution of the energy 
r 2 of the stochastic oscillator (||). Taking this formula as starting point, one 
can analytically study the transport properties of a random barrier in two 
distinct cases: the ballistic regime, when the width of the barrier is much less 
that the localisation length for the infinite lattice, and the localised regime, 
when the vice-versa is true. The two cases are respectively identified by the 
conditions L -C Zoo and loo <C L, where we use the symbol loo to denote the 



inverse of the Lyapunov exponent (16) and we are assuming that the lattice 



step is unitary, so that we can refer to L both as the number of barrier sites 
and as the length of the barrier. We will evaluate the transmission properties 
first in the ballistic and then in the localised regime. 

Before proceeding to the discussion of the two cases, we observe that our 
use of the continuous model (|5|) makes the results of this Section valid for 
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both continuous models like flI8| ) and for the discrete lattice (]]]). The same 
formulae apply to both cases, with the localisation length 1^ taking the 
form (|1|) or ( |i~6|) depending on whether the formulae refer to the continuous 
or the discrete model. We also note that the results of this section were 
obtained long ago for continuous models (see, e.g., || and references therein); 
what is new here is their application to the discrete case and the approach 
used in their derivation, which sets the mathematical results in a different 
physical perspective. 



5.1 The ballistic regime 



In the ballistic regime, i.e., when L <C l^, 
sion (E22T) can be written in the form 



one has r 12 (£) — 1 and expres- 



2 - rf(L) - r 2 2 (L) 



+ 



(23) 



Another quantity of physical interest is the resistance of the finite barrier, 
which is here defined as the inverse of the transmission coefficient 



Ri 



Tr 



2 + rf(L)+ri(L) 



(24) 



A glance at expressions ( ^3|) and ([24]) reveals that, in order to obtain the 
average value of these physical quantities, one has to compute the average 
of the squared radii rf(L) and r|(L) over different disorder realizations. To 
achieve this goal, one can rely on the method developed by Van Kampen to 
study random oscillators and other stochastic models [0, |J. Van Kampen's 
approach is based on the construction of a dynamical equation for the average 
moments of the position and momentum of the random oscillator. For the 
second moments one has 



(25) 





where the evolution matrix is 



A 





ei + e 2 
-u + e 3 





-ei + e 2 




(26) 
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with 

ei = / x(r) dr 

poo 

€2 = x( T ) cos (2o;r) dr 
e 3 = / x(r) sin (2wr) dr. 





For the general case of coloured noise, Eq. fl25|) is correct up to order 0(e) = 
0(£ 2 ); for the special case of white noise, however, it turns out to be exact. 

One can extract substantial information from Eq. (|25|); in particular, it is 
possible to obtain the behaviour of the average squared radii r\{t) and r|(t) 
for t -»• 

(r?(i)) = l + ( ei + e 2 )t + o(t 2 ) 
(r 2 2 (t)) = l + (- ei + e2 )t + o(t 2 ). 

As a consequence one has 

( r?(t)+ f ) - 2 > = ^ + ^) (27) 

Notice that this equations are correct up to order 0(t 2 ), so that it is mean- 
ingful to retain the distinction between the parameter e\ and Using the 
result ( P7|) one arrives at the following expressions for the average transmis- 
sion coefficient and resistance 

(Tl) = 1 

and 

(Rl) 

These formulae show that in the ballistic regime the averages of both the 
transmissivity and the resistance are linear functions of the thickness L of 
the disordered layer. In addition, the average resistance coincides with the 
inverse of the average transmissivity 

(TE 1 ) ^ (T L }-\ 
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5.2 The localised regime 

In the localised regime the barrier extends over several localisation lengths: 
L ^> Zoo- I n this case, to evaluate the average value of the transmission 
coefficient (|22| ) it is convenient to determine the probability distribution for 
the random variable r. We observe that for L ^> the radius increases 
exponentially and one has 

n {L) ~ r 2 (L) ~ r(L), (28) 

with probability equal to one, regardless of the initial condition. As a conse- 
quence we can drop the subscripts 1 and 2 write the transmission coefficient 
in the simplified form 

™ ^ <T7^ (29) 

From the mathematical point of view, the problem of computing the av- 
erage ( p9|) can be better handled by introducing the logarithmic variable 
z = lnr. The dynamics of the random oscillator @ is then determined by 
the equations 

* = -2^) sin2 (30) 
= uj + f (t) sin 2 0. 

System (|30D belongs to the class of stochastic differential equations of the 
form 

Ui = F/ 0) (u) + aF^ (u, t) (31) 

where (u) represents a sure function of u perturbed by a stochastic func- 
tion aFi (u, t) with a C 1. Indeed, one can reduce the system ( |3Ti| ) to the 
form fl3~T|) by defining the vectors of Eq. ([31]) as 



u 



with a — J (£ 2 (t)). It is known that a stochastic differential equation of the 



form (|31~D can be associated to a partial differential equation whose solution 



P(u, t) represents the probability distribution for the random variable u 
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This partial differential equation can be written as 



^=-i:^(^(u,*) 

+o(a 2 ) 

(32) 

where u* stands for the flow defined by the deterministic equation u = 
F^(u), d(u~ T )/d(u) is the Jacobian of the transformation u — > u~ r , and 
the symbol o (a 2 ) represents omitted terms of order higher than the second 
in the perturbative parameter a. Thus, in the case of weak stochasticity 
(q C 1), one can describe the dynamical behaviour of the system (|31"| ) with 
an approximate equation of the Fokker-Planck kind. 

In the present case, the approximate Fokker-Planck equation fl3"2|) associ- 
ated to the dynamical system (|30D takes the form 



dP dP 1 d f dP* 

-QfiO^t) = -u— + i sin(20)-|[- ei + e2 cos(20) + e3sin(2^)] — 

BP 

+ ~{ sin 2 (6) [e 3 cos(2#) - e 2 sin(20)] — 




+ sin 2 ^)^ [( ei - e 2 cos(2#) - e 3 sin(20)) P] j 
+ ~ sin(2#) [e 2 sin(2#) - e 3 cos(2#)] 

(33) 

We remark that this equation is correct to the second order in £(£) in the 
general case of coloured noise; in the special case when the noise £(t) is white, 
however, it can be shown that Eq. ( |3~3"D becomes exact. 

Once we dispose of the Fokker-Planck equation fl33|) for the general dis- 
tribution P(z, 9, t), we can consider that, in order to evaluate the average of 
the transmission coefficient (|2"2f) , we actually need only the probability dis- 
tribution for the radial variable r (or for the equivalent logarithmic variable 
z). Therefore, we do not have to solve Eq. ([33|) in all its generality and we 
can instead consider the restricted Fokker-Planck equation 

-jf P{e,z,t)d6 = -J Q [(l-cos(40))e 2 -sin(40)e 3 ] 



+ 4 



dz 2 

I r 27r dP 
/ [2 Cl cos(2#) - e 2 (1 + cos(4fl)) - e 3 sin(40) — dd 
Jo oz 



de 
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obtained by integrating Eq. ( |3"3"D over the redundant angular variable. To 
proceed further we assume that, after a short-lived transient, the probability 
distribution takes the form 

P(0,z,t)~±-P(z,t). (34) 

This assumption can be justified on the grounds that, for weak noise, the 
dynamics of the angular variable is approximately ruled by the equation 
9 ~ 10. This implies that, after a sufficiently long time (of the order of some 
periods 2tt/u), the angular variable will have swept the whole interval [0 : 2tt] 
in a almost uniform way. That makes reasonable to suppose that, for times 
t ^> 2tc/uj, the angular distribution is flat (excluding of course the exceptional 
case when u ~ 0, i.e., when the energy value lies in a neighbourhood of the 
band edge). 

As a consequence of the hypothesis (^4|), one eventually gets the reduced 
Fokker-Planck equation for the z variable 



9P , v x 



dP . s d 2 P / , 



(35) 



where A is the Lyapunov exponent fll6|) . Eq. (|35|) has the form of a heat 
equation with a constant drift; its solution is therefore 

P(z,t) = — ^exp- ^ ~ At)2 . (36) 
1 ' ' V2^Xi y 2Xt 1 ; 

This solution satisfies the initial condition P(z,t = 0) = S(z), i.e., we have 
assumed that at time t — one has r = 1, as is the case for the initial 
conditions P\ and P%. The initial condition, however, is somewhat arbitrary, 



since the equation (35) is correct only for times t ^> 2tt/lu. 

Knowledge of the distribution ( |36"D makes possible to compute the average 
transmission coefficient in the localised regime. Using probability (j3~6|) we can 



actually evaluate expression fl29|) and thus obtain 



(r '> = / r ir^y p{z - l) dz k v # ^ {- 1. 1 ■ <37) 



As a result, in the limit L — > oo one has 



~\HTl) = \. (38) 



17 



Formulae ( P7[ ) and ([JS|) show that in the localised regime the transmis- 
sion coefficient decreases exponentially with the width of the disordered bar- 
rier and they provide the correct rate of exponential decay. It must be 
pointed out, however, that expression (|37|) fails to reproduce the exact pre- 
exponential factor, which actually scales as (l tXl /L)' i l 2 (for approximation-free 
results see || and references therein). This partial shortcoming must be at- 
tributed to the two approximations made in the derivation of formula (|37D, 
i.e., i) assumption (^) that allows the substitution of the exact expres- 
sion (|22|) for the transmission coefficient with the simplified form (^) and 
ii) hypothesis (0) about the angular dependency of the probability distri- 
bution P(6,z,t). Both assumptions are admittedly incorrect for very short 
times, i.e., for distances L which are small on the length scale defined by l^. 
Thus we are led to the conclusion that in formula (|37| ) the exponential factor 
is determined by the long-time behaviour of the random oscillator (which 
is correctly described in our approach), while the pre-exponential factor is 
strongly influenced by the short-time dynamics of the oscillator. It is in- 
teresting to notice that an incorrect pre-exponential factor proportional to 
(Zoo/L) 1 / 2 was also obtained in || studying a continuous solid-state model 
with a different approach. In that study, however, the physical meaning of 
the adopted simplifying hypotheses was not so transparent as in the present 
case, where the analogy between models ([TJ) and (Q) makes possible to gain 
an intuitive comprehension of the mathematical approximations. 

Beside allowing one to compute the average of the transmission coef- 
ficient, the probability distribution (|36|) makes possible to determine the 
average value of other physical quantities which are relevant for a thorough 
description of the transport properties of a disordered barrier. The logarithm 
of the transmission coefficient and the resistance (|24D are standard choices 
for the complete analysis of the conductance problem. 

The interest for the logarithm of the transmission coefficient stems from 
the fact that, unlike the transmission coefficient itself, the logarithm lnT^ is 
a self- averaging variable and therefore a physically more sound parameter for 
the definition of the transport features of the disordered barrier (see, e.g., |]]). 
In the present framework, the average of the logarithmic transmissivity can 
be computed as follows. First, we observe again that in the localised regime 
condition ( p8|) is fulfilled for almost every realization of the disorder so that 
we can write 



18 



This expression can be put in the equivalent form 



(lnT L ) = -(In (r 2 )) + ln(2) - (In (l + 1)). (39) 

We now observe that for every x > the logarithm satisfies the relations 
< In (1 + x) < x\ hence the last term on the r.h.s of the preceding equation 
must obey 

< In (l + l)) <(!) = !, (40) 



where we have made use of distribution (pq) to evaluate the average of 1/r 



,2 



Relations ( |39"D and (fHf ) imply that in the limit L — ► oo one has 

~(^T L ) = |<lnr(L)>. 

Substituting in the r.h.s. of this equation the average value of the variable 
2 = lnr one finally obtains 

~QnT L ) = 2\ 

which shows that the average logarithm of the transmission coefficient de- 
creases linearly with the barrier width in the localised regime. 

A third quantity that represents a meaningful statistical characteristic of 
the disordered barrier is given by the inverse of the transmission coefficient, 
i.e., by the resistance fl24j) . As we did in the previous cases, we rely on the 
condition (ESI) to write the resistance in the form 

R L ^l(l + r 2 (L)). (41) 



Starting from this expression and making use of distribution (pq) we obtain 

,4L, 



w = \ 



1 + exp( 



(42) 



This expression shows that the average value of the resistance increases 
exponentially so that the resistance has a multiplicative rather that addi- 
tive behaviour as a function of the barrier length. This conclusion obviously 
ceases to be valid in the special case in which long-range correlations of the 
random potential make the localisation length diverge: in this case the 
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disordered barrier becomes transparent. We underline that, using the recipe 
given in Ref. |IJ for the Anderson model -or the prescriptions of Sec. f|for the 
continuous model (|I8D- it is possible to define a random potential such that 
the corresponding Lyapunov exponent is zero in certain frequency intervals 
and positive elsewhere. As a consequence, the disorder barrier generated 
by such a potential will be transparent for electrons with the appropriate 
energies and exponentially high otherwise. This opens the possibility of pro- 
jecting efficient electronic filters and agrees with the recent experimental 
findings discussed in 0. 

As a further consideration, we observe that Eqs. (|37| ) and fl42|) show that 
in the localised regime the inverse of the average transmission coefficient does 
not coincide with the average of the resistance 

in contrast to the ballistic regime case. 

At this point we wish to remark that the interest of expression (f|^) goes 
beyond the definition of the transport properties of a disordered barrier. 
This is so because the resistance Rl is strictly related to the energy r 2 of 
the random oscillator (|5|), as clearly shown by Eq. (pUp. The exponential 
increase of the average resistance can therefore be reinterpretated as energetic 
instability of the random oscillator (|^) on long time scales and formula ( |42| ) 
can be rewritten in the alternative form 

(r 2 (t)} oc exp (7 B i) for t > 1, 

with 

1e = 4A (43) 

where A is the Lyapunov exponent flToT) . This result shows that the energy 
of the random oscillator grows exponentially at large times (unless one has 
A = 0) and that the rate 7# of this exponential increase is equal to four 
times the Lyapunov exponent. We could have computed the energy growth 
rate also with a different approach, taking Van Kampen's equation (|25|) as 
a starting point. In fact, Eq. (^) determines the time evolution of the 
second moments of the position and momentum of the random oscillator; it 
is therefore possible to obtain the result ( |4"3"D by determining the eigenvalue 
of the evolution matrix (^) with the largest real part. 

Incidentally, we observe that some physicists working in the field of stochas- 
tic systems pretend to obtain the rate of orbit divergence (Lyapunov expo- 
nent) by erroneously dividing by a factor two the energy growth rate. The 
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mistake probably stems from (and is equivalent to) the incorrect assump- 
tion that for large times the average of the logarithm of the energy and the 
logarithm of the average energy coincide whereas the real relation is 

^<lnr 2 (i)> = ^ln(r 2 (t)> 

valid in the limit t — ► oo. 

As a last remark, we wish to point out another consequence of the corre- 
spondence between the resistance of a disordered barrier and the energy of 
the stochastic oscillator @. It is well known that in the localised regime the 
resistance Ri is a non- self- averaged quantity since the relative fluctuations 
of this quantity do not disappear in the macroscopic limit. Indeed, if we 
employ the average value fl4"2"|) of the resistance and use the distribution ( |36| ) 
to compute the average of the square of the resistance ([HP, we obtain that 
the root-mean-square deviation of the resistance behaves like 

6R L = ((R 2 L )(R L y 2 - if' 2 cc exp(2L/U, (44) 

i.e., grows exponentially with the length of the random barrier. This result 
is well known to solid-state physicists, but it may be of some interest to re- 
formulate it in terms of the dynamics of the stochastic oscillator ([5]). Then 
we can express the meaning of result ( |4*4]) by saying that the energy of the 
stochastic oscillator ([5]) is a quantity whose asymptotic value can fluctuate 
wildly from one noise realization to another noise realization. Relative fluc- 
tuations do not die for long times: in this sense it seems that the concept of 
non-self-averaging quantity can find an useful application also in the field of 
stochastic classical systems. 

6 Conclusions 

The first part of this paper is devoted to a thorough discussion of the analogies 
existing between the Anderson model ([]]) with diagonal correlated disorder 
and the stochastic oscillator (||) with frequency perturbed by a coloured noise 
with correlation function (0). Our analysis shows that the two systems are 
equivalent in the sense that there exists a close correspondence between elec- 
tronic states on one hand and space trajectories on the other. Quantitatively, 
this correspondence manifests itself in the identity of the inverse localisation 
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length for electronic states with the exponential rate of divergence of nearby 
oscillator trajectories. It is remarkable that this correspondence holds in 
spite of the fact that the Anderson model is quantum and discrete (in space) 
whereas the random oscillator is classical and continuous (in time). The anal- 
ogy between the models ([!]) and @ had already been investigated in |§ for 
the basic case of uncorrelated noise and disorder but the present work extends 
the previous conclusions taking into account the effect of correlations. 

In the second part of the work we discuss some implications of the par- 
allelism between the models (|l]) and @. In the first place, we translate 
the concept of "mobility edge" from the field of solid state physics to that of 
stochastic systems, showing how time correlations of the frequency noise may 
produce energetic stability for the random oscillator ([5]). We then show how 
knowledge about the oscillator dynamics on finite time scales can be used 
to gain insight about the transport properties of a finite disordered lattice. 
This allows us to derive significant results about electronic transmission in 
a simple and physically transparent way. Using the analogy the other way 
round, we can also deduce statistical properties of the energy of the stochas- 
tic oscillator from knowledge of the statistical features of the resistance of 
a disordered wire. In passing we also clarify the relation between energetic 
and orbit instability for the random oscillator (|5|). 

In conclusion, we believe that the bridge built among the fields of solid- 
state disordered systems and classical stochastic models represents a useful 
way to study the properties of both classes of systems. The present paper can 
be considered as an illustration of how this dual approach works, allowing 
one to solve old problems by putting them in a new perspective. 

7 Acknowledgements 

The authors are grateful to N. Makarov and V. Dossetti for fruitful discus- 
sions. L. T. would like to express his appreciation for the support offered by 
the Instituto de Fisica of Puebla, where much of this work was done. F. M. 
I. acknowledges the support by CONACyT (Mexico) Grant No. 34668-E. 

References 

[1] F. Izrailev and A. Krokhin, Phys. Rev. Lett. 82, 4062, (1999) 



22 



[2] U. Kuhl, F. M. Izrailev, A. A. Krokhin, and H.-J. Stockmann, Appl. 
Phys. Lett. 77, 633, (2000) 

[3] F. M. Izrailev, T. Kottos, and G. P. Tsironis, Phys. Rev. E 52, 3274, 
(1995) 

[4] L. Tessieri and F. M. Izrailev, Phys. Rev. E 62, 3090, (2000) 

[5] M. Kappus, F. Wegner, Z. Phys. B - Condensed Matter, 45, 15, 
(1981); B. Derrida, E. Gardner, J. Physique, 45, 1283, (1984) 

[6] I. M. Lifshitz, S. Gredeskul, and L. Pastur, Introduction to the The- 
ory of Disordered Systems (Wiley, New York, 1988) 

[7] N. G. Van Kampen, Stochastic Processes in Physics and Chemistry, 
second edition (North Holland, Amsterdam, 1992) 

[8] N. G. Van Kampen, Physica 74, 215, (1974); N. G. Van Kampen, 
Physica 74, 239, (1974) 



23 



